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I. INTRODUCTION 

The Belinskii, Khalatnikov and Lifshitz 
(BKL) scenario [1, 2] (see [3] for numerical sup¬ 
port for BKL), addresses the generic solution 
to the Einstein equations near the cosmologi¬ 
cal singularity. The purpose of our paper is to 
quantize the dynamics of the vacuum Bianchi- 
IX model that underlies the BKL scenario. 

The BKL predicts that on approach to a 
spacelike singularity the dynamics of gravita¬ 
tional field may be significantly simplified as 
time derivatives in Einstein’s equations domi¬ 
nate over spatial derivatives. The latter means 
that the evolution of the gravitational field in 
this regime is ultralocal and space splits into 
collection of small patches whose dynamics is 
approximately given by spatially homogenous 
spaces, the Bianchi models (see, e.g. [4]). Ap¬ 
proaching the singularity the spatial curvature 
grows and the space further subdivides into ho¬ 
mogenous slices. The size of each patch, mod¬ 
eled by one of the Bianchi spacetimes, corre¬ 
sponds to the magnitude of the spatial deriva¬ 
tives in the Einstein equations. As homogene¬ 
ity of spatial fragments holds only at some level 
of approximation, dynamical evolution of the 
newly formed patches starts off with slightly dif¬ 
ferent initial conditions. This almost negligible 
difference grows rapidly in subsequent evolution 
as geometries of new patches evolve almost inde¬ 
pendently of each other approaching the state of 
the so-called asymptotic silence [5] . The chaotic 
subdivisions result in the growth of fragmenta¬ 
tion of spacetime suggesting that it may be pos¬ 
sessing fractal structure close to the singularity 
[6, 7], 

Among the possible homogeneous models, the 
Bianchi-IX model has sufficient generality to de¬ 
scribe the evolution of a small patch of space to¬ 
wards the singularity. The dynamics of the vac¬ 
uum Bianchi-IX model (i.e., the mixmaster uni¬ 
verse) is nonintegrable. However, close enough 


to the singularity, each solution can be qual¬ 
itatively understood as a sequence of Kasner 
epochs, which correspond to the Kasner uni¬ 
verse. The transitions between the epochs are 
described by the vacuum Bianchi-II type evo¬ 
lution [8]. The universe undergoes an infinite 
number of chaoticlike transitions and eventually 
collapses into the singularity in a finite proper 
time [1]. 

The imposition of quantum rules into the 
chaotic dynamics of the Bianchi-IX model has 
been already studied. The program initiated 
by Misner [9-1 ] led to the pessimistic result 
that quantum mechanics does not remove the 
singularity of the Bianchi-IX model. Neverthe¬ 
less, the exploration of solutions to the corre¬ 
sponding Wheeler-DeWitt equation continues 
[12], [13], [ 4], [1 ]. Recently, some effort has 
been made towards quantization of the Bianchi- 
type models within the loop quantum cosmol¬ 
ogy. The authors make use of the Dirac quan¬ 
tization method and combine it with the in¬ 
troduction of holonomies in place of the cur¬ 
vature of connection. The results obtained for 
the Bianchi-IX model at semiclassical level by 
Bojowald [4, 16] suggest that the chaotic behav¬ 
ior stops once quantum effects become impor¬ 
tant. Another formulation taking into account 
holonomies has been proposed in [17], but it 
has not been applied to the examination of the 
dynamics. Still another proposal was given in 
[18], giving support to [16]. An effective dynam¬ 
ics considered recently [19] suggests a resolution 
of the cosmic singularity problem as well, but 
within an approach quite different from ours. 
In the above formulations, the search for solu¬ 
tions is quite challenging leaving the near big 
bang dynamics largely unexplored. 

In this paper we formulate and make a quan¬ 
tum study of Bianchi-IX model by combining 
canonical and affine coherent state (ACS) quan¬ 
tizations with a semiclassical approach. The 
cosmological system consists of isotropic vari¬ 
ables (expansion and volume) and anisotropic 
ones (distortion and shear). They are treated in 
a separate manner. The canonical pair expan¬ 
sion volume is a half plane. Since the symme¬ 
try of the latter is the affine group “AX + B ” 
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and not the Abelian M 2 , they are consistently 
quantized by resorting to one of the two uni¬ 
tary irreducible representations (UIRs) of the 
affine group. Within this approach, we have 
found in [20] that for the Friedmann-Robertson- 
Walker (FRW) models the cosmological expan¬ 
sion squared, which plays the role of kinetic 
energy of the universe, is always accompanied 
at the quantum level by an extra term in¬ 
versely proportional to the volume squared. As 
the Universe approaches the singularity, this 
term grows in dynamical significance, efficiently 
counterbalances the attraction of any matter 
and eventually halts the cosmological contrac¬ 
tion. Afterwards, the Universe rebounds and 
re-expands. In the present work (see also [21]) 
we confirm that the same mechanism prevents 
the collapse of Bianchi-IX universe, suggesting 
its universality. Making further use of the ACS 
we construct a semiclassical description of the 
isotropic part of the metric, with semiclassi¬ 
cal observables replacing the classical ones in 
the phase space. In particular, the semiclassi¬ 
cal Hamiltonian possesses the correction term, 
which regularizes the singularity. 

Inspired by standard approaches in molec¬ 
ular physics, we make an assumption about 
the quantum evolution of the anisotropic vari¬ 
ables based on the adiabatic approximation. In 
molecules, the motion of heavy nuclei is so slow 
in comparison with rapidly moving light elec¬ 
trons that it is legitimate to approximate the 
dynamics with electronic configurations being 
instantaneously and continuously adjusted to 
the position of nuclei. Analogously, we con¬ 
sider in our model, in its harmonic approxima¬ 
tion, the anisotropic oscillations rapid in com¬ 
parison with the contraction rate of the Uni¬ 
verse. Within this approach, the oscillations of 
the classical scenario are suppressed and the de¬ 
velopment of chaos is blocked. Moreover, we 
find that while the classical energy of the os¬ 
cillations behaves in terms of the scale factor a 
more or less as oc a -6 , the respective quantum 
energy behaves as oc <Y 4 , i.e. it contributes 
on a much softer level. Therefore, fluids with 
pressure equal or higher than that of radiation, 


which are likely to be present in the early uni¬ 
verse, will have their grasp on the cosmological 
collapse. 

The paper is organized as follows: Sec¬ 
tion II concerns the definition of the classical 
model, some description of its dynamics and 
the choice of phase space variables convenient 
for our quantization. Section III is devoted 
to the quantization of the Hamiltonian con¬ 
straint and its subsequent semi-classical approx¬ 
imation. Subsection IIIB explains the semi¬ 
classical Lagrangian approach. The resulting 
semiclassical dynamics is developed in Section 
IV. In Section V, we go beyond the adiabatic 
approximation to confirm the validity of our 
method. We discuss our results and conclude 
in Sec VI. In Appendix A we give an introduc¬ 
tion to the affine coherent states quantization 
together with its semi-classical aspects. Deriva¬ 
tions of the quantum version of the anisotropic 
Hamiltonian both in harmonic and triangular 
box approximations are given in Appendix B. 


II. BIANCHI-TYPE MODELS 

We consider a spacetime admitting a foliation 
M i—^ Exl, where E is spacelike. Furthermore, 
we assume E to be identified with a simply tran¬ 
sitive three-parameter group of motions. Such 
models are called Bianchi type models. The 
left-invariant vector fields are associated with 
the Killing vectors and the right-invariant ones 
with the basis vectors with respect to which the 
metric components on E take constant values in 
space. We assume the following line element 

da 2 = —Af 2 (t)dt 2 + ^ q.j (t)co l (g) a/ , (1) 

i 

where the a/’s are right-invariant dual vectors. 
They satisfy 

dw^icV^Aw*, (2) 

where are structure constants. We con¬ 
sider the so-called Class A models with C l ik = 0 
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(summation is implied). The further simplifi¬ 
cation is gained for the diagonal ones, C l j k = 
n^e^jh, where e^fc is a totally skew-symmetric 
symbol. For such models, the computation of 
the Ricci curvature is straightforward: 


where the Misner configuration variables are re¬ 
lated to the metric components as follows: 


/In 9l \ / 2 2 2^3 \ / bo \ 

lng 2 = 22 -2\/3 b+ (6) 

\ ln ?3 / \ 2 —4 0 J \j3- ) 


n 1 !! 2 rbn 3 n 2 n 3 (n 1 )” qi 

R= -+-+- 

q3 Q 2 qi 2 q 2 q 3 

(H^_ (3) 

2 qiq 3 2 q x q 2 

The vector (ni,ti 2 ,n 3 ) € R 3 specifies the 

Bianchi-type model (I, II, VIo, VIIo, VIII, IX). 
Conventionally, the rq ’s are chosen as rq £ 
{0, ±1} [22], Special cases are tq = 0 (type 
I) and rq > 0 (type IX). From now on, we fix 
= n for the Bianchi-IX case examined in this 
work. We assume that the topology of the spa¬ 
tial leaf is S 3 and we find its coordinate volume 
as: 

Vo = [ w*Aw 2 Aw 3 = (4) 

Js* n 3 

Two convenient choices are either n = 1 or n = 
\/T&7rh We make use of the latter option. 


A. Canonical formulation 


Following the work of Arnowitt, Deser and 
Misner [23], the convenient formulation of 
Bianchi models was derived by Misner [9-11]. 
With Misner’s variables the dynamics assumes 
a convenient form: motion of a particle in three- 
dimensional Minkowskian space-time and in a 
space-and-time-dependent confining potential. 
The spatial coordinates describe the anisotropic 
distortion of the shape of Universe and the time 
coordinate describes the size of Universe. The 
particle motion is ruled by a potential arising 
from the Ricci curvature of spatial leaf. Let us 
recall that the Hamiltonian constraint reads 


J\fV 0 e~ 3 ^° 

48k X 

Ho +Pl +P- - 24e 6/3o R(/3o, At)) 


(5) 


and where Po,P+ and p- are the respective mo¬ 
menta, defined from the Poisson brackets in Eq. 
7. Vo is the coordinate volume and k = 87tGc - 4 . 
The momenta pi carry the dimension L -1 while 
the positions bi are dimensionless. Because we 
reduce a field theory to a mechanical system all 
the canonical variables are in fact averaged over 
the sphere and the Poisson brackets read as: 

{bo, Po} = {b±,P±} = ^ (7) 

Fo 


In order to work within the standard quan¬ 
tum mechanical framework, we introduce vari¬ 
ables that are conjugated in the usual sense, i.e. 
{Pi,Pj} = Sij, and so define the new momenta 
Pi = (2/cc) -1 VoPi. The Hamiltonian of Eq. (5) 
now reads: 



( 8 ) 

P~) + 24e 6/3 °R(bo,b±) ) , 


In what follows we put Vo = 2k = c = 1 (in 
[21] we chose k = 1). The physical constants 
Vo, k and c will be eventually restored in the 
presentation of the final results. Note that the 
averaged (isotropic) scale factor “a” is defined 
as a = e^ 0 = {qiq 2 q 3 ) 1 ^ 6 ■ Observe also that the 
usual diffeomorphism constraints vanish identi¬ 
cally. 

In the case of Bianchi-IX geometry we find 

R(P o,b±) = -ye _2 ^+x 
^ 2cosh(2v / 3b-) - 
= -|e- 2 ^W n (b), 



( 10 ) 
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with 


W n U3) = n 2 


e 4/3 + 

~Y~ 


X 


(ii) 


^ 2cosh(2v / 3/3_) 



where the potential W n does not depend on 
the averaged scale factor a. By putting p± = 
0 = 0± we retrieve the closed Friedmann- 
Robertson-Walker model, with W„(0) = —n 2 
giving rise to the isotropic and positive intrin¬ 
sic curvature. For n = 0, we get the Bianchi 
I model with Wq = 0 and vanishing intrinsic 
curvature. 

The potential W n is bounded from below and 
reaches its minimum value, W n = —n 2 , at 0± = 
0. The potential W n is expanded around its 
minimum as follows: 

W n (0) = ~n 2 + 8n 2 (0l+0 2 _) + o(0 2 ± ), (12) 

which clearly shows its two-dimensional har¬ 
monic approximation. The W n is asymptoti¬ 
cally confining except for the following three di¬ 
rections (shown in Fig. 1), in which W n —> 0: 
(*) 0- =0, 0+ -t +oo, (ii) 0 + = 0- 

Too, (in) 0 + = ^=, 0- -T -oo. 

The form of the Bianchi-IX potential deserves 
particular attention due to its three “open” 
C 3 „ symmetry directions (see figures 1 and 2). 
One can view them as three deep “canyons”, 
increasingly narrow until their respective wall 
edges close up at the infinity whereas their re¬ 
spective bottoms tend to zero. The motion of 
the Misner particle in this potential is chaotic 
[24]. Though the curvature, which is propor¬ 
tional to the potential, flattens with time, the 
confined particle undergoes infinitely many os¬ 
cillations. In the so-called steep wall approxi¬ 
mation, the particle is locked in the triangular 
potential with its infinitely steep walls moving 
apart in time. At the quantum level, the con¬ 
fining shape originates a discrete spectrum. On 
the other hand, it is unclear on a mathemati¬ 
cal level whether or not the Bianchi-IX poten¬ 
tial also originates a continuum spectrum. Nev¬ 
ertheless the following physical reasoning leads 


to the conclusion that a continuous spectrum 
probably does not exist. Indeed a continuous 
spectrum should be associated with the eigen- 
system {(ipk(0-, 0+), £fe)}fcei> C R 2 where: (i) the 
scattering eigenstates (not R 2 -Lebesgue square 
integrable) ipk(0) represent the possibility for 
the system to go from infinity to infinity in the 
(/3_,/3_)_)-plane, (ii) k € V C R 2 represents the 
211-momentum of the incoming “particles”, (iii) 
k belongs to some 2Z)-domain V of R 2 (T> is not 
restricted to a subset of measure 0 for the R 2 - 
Lebesgue measure), (iv) k G V leads to an inter¬ 
val (generally unbounded) of eigenvalues (the 
continuous spectrum). In our case the potential 
only possesses three discrete “open” directions 
at infinity, therefore the incoming momentum k 
cannot belong to a 2D-domain V (a subset with 
a non-vanishing R 2 -Lebesgue measure). Then 
it is reasonable to think that a continuous spec¬ 
trum does not exist. 



FIG. 1. Contour plot of W„(0) near its minimum. 
Red lines: the three C 3 „ symmetry axes 0- = 0, 


0+ = 0./V3, 0+ = -0-/V3. 


The evolution of Bianchi IX can be viewed as 
a non-linear model of gravitational wave in dy¬ 
namical isotropic geometry (see e.g. [10, 25] and 
references therein). The wave, which consists of 
two nonlinearly coupled components 0±, is ho- 
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mogenous, that is, it is a pure time oscillation. 
Its wavelength is thought to be much larger than 
the extension of the considered patch of the 
universe with its spatial derivatives neglected. 
The energy of the wave sources the gravitational 
contraction. As we later show, the quantiza¬ 
tion of the wave introduces important modifi¬ 
cations to the dynamics of the whole universe. 
The qualitative study of spatially homogeneous 
models was pioneered by Bogoyavlensky (see 
[26] and references therein). 


B. Redefinition of basic variables 


For the purpose of ACS quantization, we re¬ 
define the isotropic phase space variables as the 
canonical pair: 

q '■= e 3/3 °/ 2 , p:=\e~ 3 M 2 p c (13) 

Note that ( q,p ) lives in the half-plane. The 
Hamiltonian (8) now assumes the form 


H = - 


24 




36q 2/3 W n (p)^J . 


(14) 

Let us split the potential W n into its isotropic 
and anisotropic components. The isotropic part 
corresponds to that part of curvature which is 
independent of /3 ±, whereas the anisotropic part 
vanishes for /3± = 0: 


Wn(P) = -n 2 + V n (P) (15) 


The lapse function is not dynamical and its 
choice is irrelevant for the classical dynamics 
as it only fixes the magnitude and direction of 
the Hamiltonian flow in the constraint surface. 
From now on we put the lapse A f = —24. The 
Hamiltonian (14) reads now 

H = jjp 2 + 36nV /3 -H 9 , (16) 

where H g is the g-dependent Hamiltonian for 
the anisotropic variables, 

y^ + 36 <? 2/3 K03). (17) 


C. Discussion of the constraint 

The analytical expressions (16)-(17) for H 
remind us of molecular system’s Hamiltonian. 
The pair (q,p) plays the role of the nucleus-like 
dynamical variables and (j3± , p± ) are electron¬ 
like dynamical variables. Quantum molecular 
systems are usually considered by making use 
of the Born-Oppenheimer approximation (BO) 
or its “diagonal correction” named Born-Huang 
(BH) [27, 28]. In molecular physics, the valid¬ 
ity of these approximations depends crucially 
on the ratio between nuclei and electron masses. 
Namely, a nucleus mass is very large when com¬ 
pared to the electron mass. In our case, Eq. (17) 
indicates that q 2 plays a role of “mass” for the 
anisotropic variables, whereas Eq. (16) shows 
that the “mass” of the isotropic variable is con¬ 
stant. Therefore, near the singularity q = 0, 
P± ’s become light degrees of freedom and q with 
its constant mass may be treated as a heavy 
degree of freedom. Hence, we may follow ei¬ 
ther the BO or the BH approximation scheme 
in quantizing our system. This issue is consid¬ 
ered in more detail in the next section. 

For Eqs. (16)-(17), one checks that /3± = 0 = 
p± is a solution to the Hamilton equations of 
motion. In this case, the constraint H = 0 re¬ 
duces to 

^p 2 + 36n 2 g 2/3 = 0. (18) 

and we recover the closed vacuum FRW con¬ 
straint, which possesses the unique singular and 
uninteresting solution p = 0 = q. Nevertheless, 
it makes sense to consider a small perturbation 
SP± from i3± = 0. The dynamical equation for 
S8± based on the harmonic approximation of 
Eq. (12)) is 

. 2 

= -2 - 2 . (19) 

q g 4 / d 

It can be demonstrated that the Friedmann 
model evolving towards the singularity is not 
stable in the phase space of the Bianchi IX 
model. More precisely, any perturbation of 
isotropy will grow and develop into an oscil¬ 
latory and chaotic behavior. The growth of 
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isotropy is apparent from Eq. (19) for q/q < 0. 
As the shear grows, the harmonic approxima¬ 
tion breaks down and fully non-linear dynam¬ 
ics develops. As we show, this behavior is sup¬ 
pressed on the quantum level and allows the 
harmonic approximation to remain valid all the 
way towards the big bounce. Moreover, we 
show in Section V that the BO approximation 
also survives the bounce phase. As the shear 
is known to dominate over any type of famil¬ 
iar matter close enough to the singularity, we 
consider only the vacuum case. Then the lack 
of solutions for the vacuum isotropic universe 
as concluded from (18) leads to the prediction 
that one should take into account the effect of 
the ‘quantum zero point energy’ generated by 
the quantized anisotropy degrees of freedom of 
the Bianchi-IX model. 


V,0S_,fi + ) 



p. 

FIG. 2. The plot of V„ for n = 1 near its mini¬ 
mum. Boundedness from below, confining aspects, 
and C 3v symmetry are illustrated. 


III. QUANTIZATION 

In what follows we apply a quantization 
based on the Dirac method and inspired by 
Klauder’s approach [29]: (i) quantizing H 


in kinematical phase space, (ii) finding the 
semi-classical expression H of the quantum 
Hamiltonian H using Klauder’s approach and 
our adiabatic approximation, and (iii) imple¬ 
menting the Hamiltonian constraint on the 
semi-classical level H = 0. 


A. Quantum constraint 

Since in Eqs. (16)-(17) we have ( q,p ) G R+ x 
R and (/ 3±,p ±) £ R 4 , we can follow the idea of 
our previous paper [20] in the realization of the 
step (i): 

(i) for the quantization of functions (or dis¬ 
tributions) of the pair ( q,p ) living in the 
open half-plane, we apply ACS quanti¬ 
zation, whose principles and methods, as 
part of integral quantizations (see [30] and 
references therein), are explained in [20] 
and summarized in Appendix A. This pro¬ 
cedure yields p = — ihd x and q defined as 
the multiplication by x , both acting in the 
Hilbert space L 2 (R + , dx ), 

(ii) for quantizing functions of the pairs 
{(3±,p±) we can choose either the Weyl- 
Heisenberg coherent state (CS) quantiza¬ 
tion [30], which has a regularizing role, 
or, directly, the canonical quantization. 
Due to the simplicity of the model we fol¬ 
low the later option. Actually both yield 
p± = — ihdp ± and the multiplication oper¬ 
ator /3 ± , both acting in L 2 (R 2 ,d^ + d/3_). 

Thus, for the quantum Hamiltonian H corre¬ 
sponding to Eqs. (16)-(17) we get 

H = \ if + ~<^) + 36n2 ^ 2/3 ~ H? > ( 2 °) 

H, := # /+ \ P - + 36£ 3g 2 / 3 U n (/3) ■ (21) 

q- 

where the are positive and non-vanishing 
purely numerical constants dependent on the 








B. Semi-classical Lagrangian approach 


choice of the so-called fiducial vector. With 
the choice made in our previous paper [20] , and 
thanks to the formula recalled in Appendix A, 
we have 


•&1 —4 ( 1 + v 

£3 = 


Kq{v) 
Ki{v). 

-^5/3 M 


^2 = 


K 2 (y) 

Ki(v) 


K 1 (uy/ 3 K 2 { v)V 3 


( 22 ) 


where v > 0 is a free parameter and the K r {y) 
are the modified Bessel functions [32]. Since we 
deal with ratios of such functions throughout 
the sequel, we adopt the convenient notation 


frs — frs (Z'') 


K r {y) 

K s (y) 


1 

f,sr 


(23) 


One convenient feature of such a notation is that 
frs{v) ~ 1 as v —> 00 (a consequence of K r {y) ~ 
y/wj (2 1 /)). Thus (22) reads 


£i = | (i + i'&iM) , ^2 = (6i(^)) 2 , 

/ \!/3 / \ 2/3 

^ = (^fi) (^fa) ( 24 ) 


There exist many other choices of fiducial vec¬ 
tors yielding similar constants A.;. These vec¬ 
tors depend themselves on arbitrary parame¬ 
ters which can be suitably adjusted. Also, the 
Hamiltonian (20) itself is defined up to a mul¬ 
tiplicative factor. It is also crucial to recall 
that there exists an infinite range of values for 
constants A,;, for which the symmetric H has 
a unique self-adjoint extension. This is proved 
by making use of the reasoning previously pre¬ 
sented in [20] and recalled in Appendix A. 

Important remark: The analysis developed in 
the following sections is independent of the spe¬ 
cific numerical values given to the coefficients 
as long as the bounce happens for small 
enough values of q in order for the adiabatic 
approximation to be valid (see below). In other 
words, all the following results determine the 
qualitative features of the dynamics unambigu¬ 
ously and the numerical coefficients pin down 
the dynamics quantitatively. 


Being inspired by Klauder’s approach [29], we 
present a consistent framework allowing to ap¬ 
proximate the quantum Hamiltonian and its as¬ 
sociated dynamics (in the constraint surface) by 
making use of the semi-classical Lagrangian ap¬ 
proach. 


1. General setting 

The quantum Hamiltonian (20) has the gen¬ 
eral form (up to constant factors) 

H =M[f+ I ^+Lqi , (25) 

where K and L are some positive constants and 
the (/-dependent Hamiltonian H( int )(</) (also de¬ 
noted by H (; in the previous subsection) acts on 
Hilbert space of states for ‘internal’ degrees of 
freedom, i.e., the anisotropic ones. 

The Schrddinger equation, ih^\^{f)) = 

A/’H|H/(t)), can be deduced from the Lagrangian: 

L(*,vi/,A0 :=<*(t)| !*(*)>, 

(26) 

via the variational principle with respect to 
|T(t)). In order to solve the Schrodinger equa¬ 
tion, the lapse function A f should be fixed. 

The quantum counterpart of the classical con¬ 
straint H = 0 can be obtained as follows 

= {H,(t)\p 2 +^+Lqi-H^(q)Mt)) =0. 

,(27) 

The commonly used Dirac’s way of imposing a 
constraint, H| v k(t)) = 0, implies (27) but the 
reciprocal does not hold in general. 

At this stage, we suppose (due to the con¬ 
fining character of the potential V n ) that there 
exists H( int ) (g) as self-adjoint operator (and as 
a function of the c-number q) with purely dis¬ 
crete spectral decomposition, which is of course 
true in its harmonic approximation, 

H (int) ( 5 ) = £ E^\q) \<t>m<t>n • (28) 
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To present the Klauder semi-classical pro¬ 
cedure in the most general situation (not re¬ 
stricted to Bianchi-IX), we distinguish the two 
cases: 

(i) is independent of q , which allows 
a complete separation of variables, and 
leads to the original Born-Oppenheimer 
[27, 28] approach; 

(ii) is dependent on q. 

Different semi-classical approximations result 
depending on these choices 1 * . 

2. Semi-classical Lagrangian approximations 

(i) ^n int) independent of q 

In this case, a family of exact solutions 
of the time-dependent Schrodinger equation 
= H|(t)) can be introduced in the 
form of the tensor product 

!*(*)> = \Ht)) <8 | 4 int) >, ( 29 ) 

where | <fi(t)) is solution to the reduced time- 
dependent Schrodinger equation 

*^10(0) = ( 3 °) 
a f (p 2 + ^_+ - ^ int) (9)) \m =■■ 

where En nt is the eigenvalue of H^"). The 
tensor product |\I/(f)) of Eq. (29) is precisely 
the Born-Oppenheimer-like solution [33]. The 
equation (30) may be derived from a variational 
principle applied to the quantum Lagrangian 

L lod (<^,A0 := (31) 

i m- 


1 We first present the case (i), which is simple, and later 

introduce the more complicated case (ii), being ap¬ 

plied to the Bianchi IX model. 


Following Klauder [29], we assume that \cj)(t)} 
is in fact an affine coherent state. We as¬ 
sume in the following that the fiducial vec¬ 
tor if) has been chosen such that Co(ip) = 
in order to obtain the canonical rule 
[A g ,A p ] = ih. Furthermore, we need to apply 
a rescaling \q{t),p(t)} —► \Xq(t),p(t)} in order 
to ensure (Xq{t),p(t)\ A q \q{t),p(t)} = q(t) and 
(Xq(t),p(t)\A p I q(t),p{t)) = p{t). The parame¬ 
ter A is uniquely defined by the choice of the 
fiducial vector, namely A = 1 /c _3 (V 7 ) ( see ap¬ 
pendix A). 

Therefore we replace [/['(t)} in L of Eq. (26) 
by 

\*it)) = \Xq{t),pit))®\^). (32) 

where q(t ) and p(t) are some time-dependent 
functions. Then the Lagrangian (26) or (31) 
turns to assume the semi-classical form 

L sc (q,q,P,P,Af) = (33) 

{Xq{t),p(t)\ (ih^ - A/’H( t ed ^ \Xq(t),p{t)) 

= -qp-Af(Xq{t),p(t)\f\™ d \Xq(t),p(t)} 

= +dP 

at 

—A/" {Xq(t),p(t)\\-\™ d \Xq(t),p{t)). 

The appearance of the first term — qp in the 
r.h.s. of this equation results from the deriva¬ 
tive of (A12) with respect to parameters q and 
p leading to (A31). 

The semi-classical expression for the Hamil¬ 
tonian is the lower symbol 

H™ d (q,p) := (Xq,p\H™ d \Xq,p) . (34) 

It is defined by the ‘frozen’ quantum eigenstate 
‘n’ of the internal degrees of freedom. 

From this reduced Hamiltonian one derives 
the equations of motion together with the con¬ 
straint 


q = A/'^H(T d (g,p), 

(35) 


(36) 

0 = H r : d (q,p). 

(37) 
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These equations will allow us to set up 
Friedmann-like equations with quantum correc¬ 
tions for g and p. 

(ii) dependent of q 

Let us examine the general case in which the 
eigenstates |</>n nt ^) depend on q. We start again 
from the spectral decomposition 

H (int) (<z) = J2 Ei n nt) (Q) \^ nt) (Q)M nt) (Q)\, 

n 

(38) 

and we pick some other g-independent orthonor¬ 
mal basis |e£ nt) > of the ‘internal’ Hilbert space 
. This change of basis is associated with 
the introduction of the g-dependent unitary op¬ 
erator 

(39) 

n 

which allows to deal with the analogue of the 
Hamiltonian (28): 

H (int) (g) =U\q)H (iat \q)U(q) = (40) 

n 

The quantum Hamiltonian (25) has now the 
general form 

H = Af (V + | + Lgl - U(q)H^\q)U\q)^j . 

(41) 

The difference between the Hamiltonians 
of cases (i) and (ii) is the presence in (ii) of 
the unitary operator 17(g) that introduces a 
quantum correlation (entanglement) between 
the ‘internal’ degrees of freedom (anisotropy) 
and the ‘external’ one (the scale factor). 
As a consequence, any solution |'L(f)) of 
the time-dependent Schrddinger equation 
iftj^l'l'(f)) = H|’F(t)) cannot be factorized as a 
tensor product like ]</>(£)}<g>|<^ lnt 9(t)), contrarily 
to the case (i). In our case we wish to follow 
Klauder’s approach to build some semi-classical 
Lagrangian analoguous to Eq. (31). We use 
the previous case (i) as a starting point (as a 


guide) to build approximate possible forms of 

I *(*)>• 

It is interesting to notice that the Hamil¬ 
tonian H of Eq. (41) is unitarily equivalent to 
the one that occurs in quantum electrodynam¬ 
ics. For this purpose let us introduce the q- 
dependent operator A(g) acting on the Hilbert 
space of internal degrees of freedom as 

Mq)=ih^( q )U\q). (42) 

As a matter of fact, A(g) is self-adjoint, and the 
Hamiltonian H of Eq. (41) reads as 

H = AfU(q)x (43) 

((P - A(g)) 2 + p + Lqi - H< int )(g)) U\q). 

The form (43) (i.e. modulo the unitary trans¬ 
formation 17(g)) is the most general form of 
the Hamiltonian of a particle moving on a 
half-line and minimally coupled to an external 
field. This form is also most general on the 
grounds of the so called shadow principle [34]. 
According to the principle, at any fixed time 
the state of an interacting particle should 
be indistinguishable from the state of a free 
particle and the Galilean addition rule of 
velocities should be preserved. In particular, 
if we interpret H( mt )(g) as a g-dependent elec¬ 
tromagnetic energy (despite all the differences) 
and A(g) as a gauge field, the problem appears 
(up to a unitary transformation) similar to the 
one of a charged particle in interaction with an 
electromagnetic field. 

Now, using (41) and taking into account the 
analysis of the previous case (i), one can define 
different possible expressions of |4/(t)): 

(a) In the first approach we keep the tensor 
product expression of (i), but inserting 
the g-dependence of eigenstates. This cor¬ 
responds to a Born-Oppenheimer-like ap¬ 
proximation: 

\^{t))^\\q{t),p{t))®\^\q{t)))- ( 44 ) 
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(b) The second strategy consists in introduc¬ 
ing some (minimal) entanglement between 
q and ‘internal’ (anisotropy) degrees of 
freedom. This corresponds to a Born- 
Huang-like approximation: 

|^(t))«C/(g)(|Ag(t),Ki))®|el int) )) . (45) 

(c) In the third method one keeps the ten¬ 
sor product approximation, but including 
a general time-dependent state for the in¬ 
ternal degrees of freedom: 

|T(t))«|Ag(t),p(t))®|^(t)). (46) 

(d) The fourth strategy is the most general 
one. It consists in merging (b) and (c): 

|vk(t))«t/(g)(|Ag(t),p(f))®|0( int )(t))) . 

(47) 

Building now the semi-classical Lagrangian in 
agreement with the procedure defined in (i), we 
can distinguish two categories in the approxima¬ 
tions listed above. 

(1) (a) and (b) are completely manageable 
on the semi-classical level: they involve 
q and p as dynamical variables, while the 
anisotropy degrees of freedom are ‘frozen’ 
in some eigenstate; (a) and (b) correspond 
typically to adiabatic approximations. 

(2) (c) and (d) are more complicated: 
they mix a semi-classical dynamics for 
(q,p) and a quantum dynamics for the 
anisotropy degrees of freedom. This 
corresponds to ‘vibronic-like’ approxima¬ 
tions, well-known in molecular physics 
and quantum chemistry [35]. In our case 
this means that different quantum eigen¬ 
states of the anisotropy degrees of freedom 
are involved in the dynamics: during the 
evolution excitations and decays are pos¬ 
sible, with an exchange of energy with the 
‘classical degree of freedom’ (q,p). 


The Bianchi-IX Hamiltonian belongs to the 
general case (ii). So the different approxima¬ 
tions (a), (b), (c), (d), presented above can be 
tested. In this paper we restrict ourselves to the 
presentation of the simplest cases (a) and (b). 
We postpone the study of the more complicated 
cases (c) and (d) to future papers. 

C. Semiclassical constraint 

We now proceed to the spectral analysis of 
the operator H by making use of BO-like and 
BH-like approximations, presented in a general 
form in subsection III B, and to its semi-classical 
analysis through affine coherent states. 

1. Born-Oppenheimer approximation 

In this approximation we assume that the 
anisotropy degrees of freedom are frozen in some 
eigenstate | (f>n nt \q{t))), evolving adiabatically, 
of the g-dependent Hamiltonian H ? given by 
(21). If we denote by E^{q) the eigenenergies 
of Hq, the reduced Hamiltonian H^ d of Eq. (30) 
reads 

H N d = | [f + + 36n 2 H 3 g 2/3 - E N (q). 

( 48 ) 

Due to the harmonic behavior of V n near its 
minimum, i.e., 

W) = 8n 2 (Pl+P 2 _) + o(Pl), (49) 

the harmonic approximation to the eigenener¬ 
gies Eu(q) is manageable (N = 0,1,...), giving 

E N (q) ^ ^n^/2& 2 & 3 (N + 1). (50) 

In fact N = n + + n_, n± G N. More details, 
including a discussion of another common ap¬ 
proximation made for the anisotropy potential, 
are given in Appendix B. 

Remark: The harmonic form of E^{q) is an 
increasingly rough approximation for large 
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values of N, since V n is highly non-harmonic 
far from its minimum. But for small values of 
N, this expression is valid for any value of q. 
The steep wall approximation (see Appendix 
B) is able to give a better expression for the 
eigenenergies as their values go to infinity. 
However, both approximations do not change 
the main line of reasoning in what follows. 

Taking into account the rescaling of affine 
coherent states (see Sec. Ill B 2) , the semi- 
classical expression H^ d involved in Eq. (34) is 
defined as 

H^ d (g,p) = (Ag,p|H^ d |A g ,p), (51) 

where A := £ 0 2 (v) is chosen to get the exact 
correspondence (see Sec. Ill B 2) 

(\q,p\q\\q,p) = q, (Xq, p\p\Xq,p) = p . (52) 

Finally, we obtain 

H N d {<hP) = | (p 2 + h ~^j + ( 53 ) 

36n 2 A 5 g 2/3 - + 1), 

where the three new constants A, = A(z') are 
given by 


anisotropy degrees of freedom are frozen in some 
eigenstate, but it is an eigenstate |en Ilt ' ) ) of the 
Hamiltonian bh lnt l(g) introduced in Eq. (40). In 
the case of Bianchi-IX, thanks to the harmonic 
approximation of V n (f3), we can find an approx¬ 
imation for H ,lnt \q) and the corresponding uni¬ 
tary operator U(q) in Eq. (39). We get 

17(g)~eTr (ln9)f) , (56) 

with 

D = D+ + D_, D± = t^(p±/3± + P±P±) ■ 

(57) 

We deduce for the ‘gauge field’ A (q), the ex¬ 
pression in Eq. (42), becomes 

k(q) = - 2 ^ D. (58) 

3 q 

As it is shown in the appendix B, the har¬ 
monic approximation implies that the eigen¬ 
states |ei mt) ) depend in fact on two positive 
integers n± (we use the notation lei 1 "^)). It 
results 

(el“ t) |D|el“ t) )=0, (59) 

(e(“)|D 2 |e( i f)) = 2 « + n 2 _ + n + + n_ + 3). 


H4 = (£io) 2 (jj; + ^61 + -g-£io + ^ 20 ^) > 

/ \ 1/3 / \ 2/3 

£ 6 =\/2 ( Ci 0 ) 4/3 ( 6 o ) 5/3 ( C | i ) 1/2 ■ 

(54) 

For large values of v (typically v > 20) we get 
A 4 -^, As - 1, A 6 -y2. (55) 


Therefore, the expectation value ('k(t)|H|4'(t)) 
of the Hamiltonian H, as defined in (43), for 
j'F(f)) = T r (g)(|Ag(t),p(t))(g)|ei I " t) )) correspond¬ 
ing to the case III C 2 (Born-Huang-like approx¬ 
imation), reads 


<*(t)|H|*(t)> = 

(Xq(t),p(t)\ ^ 


(60) 


-2 + x(n+,n_)) 

P ' ±0 


36n 2 A 3 g 2/3 - E N (q) \Xq(t),p(t)) 


2. Born-Huang approximation 

In the Born-Huang approximation framework 
(see Sec. IllB 2), we also assume that the 


with 

N=n + +n—, and (61) 

2 

X(n + ,n_) = -{n\ +n 2 _ + n+ + n_ +3). 
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Hence, up to the modification i —> + 

X(n + ,n_), we recover the previous expression 
of Eq. (48) for the Born-Oppenheimer-like case 
IIIC 1. From Eq. (53) we deduce the final ex¬ 
pression of the senri-classical Hamiltonian (in 
the harmonic approximation of V n (J3)) 

H r n e ± d fep) = (62) 

9 / 2 . S 2 (J?4 + i^x(n+,n_))\ 

+ -f-; + 

36n 2 ^ 5 g 2/3 - + 1), 

where £4 = (£10 f £20 ■ 


IV. SEMICLASSICAL DYNAMICS 

The quantum-corrected constraint 
= 0 may be interpreted as a 

semiclassical version of the Friedmann equa¬ 
tion. The anisotropic degrees of freedom, 
averaged at the quantum level, give rise to the 
isotropic radiation energy density. This energy 
gravitates like common matter and fuels the 
contraction. It leads to a supplementary term 
in the Friedmann equation, shown below. 


A. Effective Friedmann-like equation 


Rewritten in terms of the scale factor a, the 
constraint H a v{<hP) = 0 reads 


— , ~ 2 /, 2 1^4 87 tG . . 

— + k— + c 2 l 2 — = —p(a), 


a 2 a 2 a 6 

where, using the Planck area ap = 2ttGHc~ 3 , 
•S.5n 2 „ V, 


(63) 


k =- 


'0 
a p' 


(64) 


p(a) = nV 0 


hc(N + 1) 


The classical constraint (18) is recovered for 
h —> 0. The main features of this quantum cor¬ 
rected model are: 


(i) The value of the isotropic curvature, 
kc 2 a~ 2 , present in closed FRW models, is 
dressed by the quantization with a con¬ 
stant £ 5 , which is close enough to 1 to be 
ignored in qualitative considerations. 

(ii) The repulsive potential term proportional 
to a~ 6 , absent in classical FRW/BIX 
models, is generated by the affine CS 
quantization. 

(iii) The energy of the anisotropic oscillations 
is turned at the quantum level into the 
radiation energy, p{a). The expression for 
p(a) in terms of the quantum number N 
becomes a poor approximation for high 
values of N , due to the breakdown of har¬ 
monic approximation. Nevertheless, the 
dimensional analysis shows that the de¬ 
pendence p(a) oc nV^ 1 hca~ A is correct 
for a —> 0 . 


B. Comparing classical with semiclassical 
constraint 


The classical Friedmann equation generalized 
to Bianchi IX geometry reads 

%+ l3 R- 1 S 2 = 0, (65) 

a- o 0 

where 3 R and E 2 are respectively the curvature 
and shear of the spatial sections. The Bianchi 
IX curvature (9) may be splitted into isotropic 
and anisotropic parts and the shear may be ex¬ 
pressed in Misner’s canonical variables as fol¬ 
lows: 


3 R = 3 


1 - V{P) 

2 a 2 ’ 


,2 = P+ +P- 
24 a 6 


( 66 ) 


where the anisotropic part of curvature poten¬ 
tial V(/3), defined by ( 11 ) and (15), reads ex¬ 
plicitly as: 


e 4/3+ 

n/3) = — x (67) 

( (V 6/3 + - 2 cosh(2v / 3^_)) 2 - 4J + 1. 
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This way we arrive at the following form of (65) : 


a 2 1 _ 1 P+ +P- V(0) 
a 2 4 a 2 6 24a 6 4a 2 


( 68 ) 


The left-hand side of (68) contains only the 
isotropic variables, whereas the right-hand side 
contains all the anisotropic terms. The first 
and the second term of (63), namely a 2 a -2 and 
kc 2 a -2 , correspond to the first and the second 
term of (68), i.e. to a 2 a~ 2 and (2a) -2 , respec¬ 
tively. The third term in (63), c 2 £ -2 J?4a -6 , has 
no corresponding term in the classical equation. 
It is a purely quantum effect induced by the 
affine CS quantization. The right-hand side of 
(63) describes the quantized anisotropic energy, 
which depends on the scale factor as a -4 . It 
corresponds to the energy of anisotropy oscil¬ 
lations in the r.h.s. of (68), which depends on 
the scale factor roughly as a -6 . This difference 
between the right-hand sides of (63) and (68) 
in dependence on a is another purely quantum 
effect. 

Concerning the curvature-like term 
H5n 2 c 2 (2a) -2 in the semiclassical Eq. (63), we 
note that it reproduces the classical re-collapse 
in the evolution of the vacuum Bianchi IX 
model, which results from S 3 topology of this 
model [31]. The re-collapse is included in the 
plots of Fig. (3). The fact that our semiclassical 
procedure reproduces this particular classical 
feature can be viewed as the confirmation of 
the consistency of our approach. 

Furthermore, from the semi-classical La- 
grangian approach developed in Subsection 
IIIB, we note that any effective matter term 
F(q) = p(q 2 / 3 ) with p(a) = Ka~ n may be in¬ 
corporated in the Lagrangian of Eq. (25), lead¬ 
ing after the action minimization of Eq. (27) to 
a supplementary matter term p(a) = £ Ka~ n = 
£p(a), (£ > 0), in the r.h.s. of Eq. (63). There¬ 
fore adding an effective radiative matter term 
like po° -4 (where po is a positive constant) into 
the r.h.s. of (63) would bring only a quantita¬ 
tive change in the semiclassical dynamics as it 
scales exactly as quantized anisotropy. 


C. Singularity resolution 


Equation (63) implies 

kc 2 V c 2 £-' 2 '^ - ^Y~a 2 p(a) < 0, (69) 

CL OC 

which defines the allowed values assumed by a. 
The inequality (69) can be satisfied only for 

< w+1 > 2 ^Tp < 70 > 

or (N + l) 2 > /(i/) with 

f( v ) := 32 (^ 02 ) 4 (£ 10 ) 1/3 £| lX ( 71 ) 

(v 2 v c 3p t \ 

Vl6 + 4 e21 + ¥ a ° + 6 °J ' 

The function f(v) is strictly increasing, from 0 
to Too, with f(u) — 1/2 i n the limit v —* oo. 

Therefore for each value of N, there exists a 
bounded interval ]0, v m (N)] for which the condi¬ 
tion (71) holds true, with = (N + l) 2 . 

For large values of N, we have v m {N) — 
^(iV+1), and u m (N) > i/ m (0) ~ 7.4. 

Thus, for v €]0,i/ m (0)[ and for all N > 0, 
the condition (69) is satisfied. We find that 
a £ [a_, a+], where 


8ap Hg 

3nV 0 it 5 


(N+l) 


1 ± 



fH 

{N Tl) 2 


(72) 

Therefore, the semiclassical trajectories turn 
out to be periodic, and a is bounded from be¬ 
low by a_. This demonstrates that the system 
does not have the singularity that occurs at the 
classical level of the FRW/BIX model. We note 
that a± is dimensionless as nVo <x n -2 is ho¬ 
mogeneous to an area. Moreover, the volume 
of the universe, a 3 Vo, is independent of n. This 
proves that our result is physical. 

In Fig. 3 we plot a few trajectories in the half¬ 
plane ( a,H). The classical closed FRW model 
is recovered at H = 0 and large values of v. 

Some remarks: 
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FIG. 3. Three periodic semiclassical trajectories 
in the half-plane ( a,H ) from Eq. (63). They are 
smooth plane curves. We use standard units ap = 
c = h = 1 and choose v = 3, Vo = 1 (so 
n = (167T 2 ) 1 / 3 ). Blue dotted curve for N = 0, green 
dotdashed for N = 1 and red dashed for N = 2. 
Each plot includes the quantum bounce (at small 
a) and classical recollapse (at large a). 


(i) The product a_a+ is only dependent on 


2a_p / J ?4 

+ r>Vo V 


(73) 


(ii) Specifying TV and v = v m (N), the relation 
a_ = a+ holds true so we have an unusual 
feature of a stationary universe with finite 
radius. 

(iii) For TV = 0 and v e]0, ^ m (0)[, the model 
shows the effect of the ‘quantum zero 
point energy’ of the anisotropy degrees of 
freedom. 


where E is the complete elliptic integral 
of the second kind [3 ]. 


V. BEYOND ADIABATIC 
APPROXIMATION 

As we have seen in Section (II), when the 
isotropy of the closed FRW universe is per¬ 
turbed, the universe acquires the Bianchi IX ge¬ 
ometry and a small perturbation inevitably de¬ 
velops into the chaotic regime, first described in 
Ref. [1], To see whether such behavior is possi¬ 
ble at the quantum level, we now go beyond the 
adiabatic approximation and allow the number 
of particles TV to grow as the universe contracts, 
bounces and re-expands. We assume that the 
scale factor is a c-number, which evolves ac¬ 
cording to the semiclassical constraint (63) for 
TV = 0. The anisotropy degrees of freedom are 
quantized as before but with the possibility to 
be excited by the time-dependent background. 

We resort to the well-known fact about the 
harmonic oscillator that its classical and quan¬ 
tum dynamics of the basic variables are in one- 
to-one correspondence and we will work in the 
Heisenberg picture by solving classical equa¬ 
tions of motion. First, we find the semiclassi¬ 
cal evolution of a in suitable time parameter, r. 
Then we use this evolution to define our wave in 
the time-dependent background in terms of an 
oscillator with “time-dependent mass”. Finally, 
we perform numerical computations. 


A. Evolution of scale factor 

The semiclassical constraint (63) may be 
written as: 


a 6 a 2 a 4 ; 


(75) 


(iv) For v e]0, v m { 0)[ and all TV > 0, the oscil¬ 
lation period T of the universe is 


T = 


4 

ncv^Ss 


a_ E 




(74) 


In what follows we set Ao = 0. This simplifica¬ 
tion removes the classical re-collapse, the uni¬ 
verse expands forever now. This is a good ap¬ 
proximation as we are not interested in the clas¬ 
sical phase but in the quantum one during which 






16 


the isotropic intrinsic curvature is assumed neg¬ 
ligible. Moreover, we may work in natural units 
in which Ai ~ A 2 as both terms origin from 
quantum theory. 

We find from (75) 



where t is the cosmological time. 

The dynamics of anisotropy is given by 
background-dependent Hamiltonian: 

H± = \p\ + 144nV(r)/3 2 ( 77 ) 


which is a part of the Hamiltonian constraint 
(16)-(17) if the lapse is set to M = —12a 3 . The 
idea of the subsequent calculations is to treat 
a(r) as a fixed function of time. To obtain the 
correct solution a(r), we need to adjust (76) for 
the choice of lapse M. We find 


dr = 


d t 

Jf 



hence 


12VA^ 


arcsm 



a(r) = 


sin ( 12 VA 2 t + |) 


(78) 


(79) 


where r € (,-247W’247 wJ &nd ° G 



B. Excitation of quantum oscillator 

The Hamiltonian under study is 

H± = + ^w 2 (r)/% (80) 


reads: 


^ = -w 2 (t)/3 (81) 

We will work in the Heisenberg picture, and as¬ 
sume that 

P{r) = {av* (t) + a) v{t)) , (82) 

p(t) = ^ (av*{r) + o*v{t)) 

where a and a ' are fixed operators and where 
v(t) solves the equation (81), i.e., 

^ = —uj 2 (t)v (83) 

We demand the canonical commutation relation 
and we obtain 


H = \P,p\ = ( av*(r) + a t n(r)) , (84) 

(av*(r) +a f v(r))] = [a, a 1 ’] V V ? VV 

We find from e.o.rn. (83) 

(v*v — vv*) = 0 (85) 

dr 

and we fix v*v — vv* = 2 i. So, a and a) are 
annihilation and creation time-independent op¬ 
erators. All time dependence lies in v(t). 

The Hamiltonian now reads 

H := \f + \^T)P (86) 

= £ ««*f+6/V) 2 ) + + „ V ) + 

We set the vacuum state |0) for a, a' to minimize 
the energy at some initial moment ro- It follows 
that 


where u> = 12y2na 2 (r). In what follows we 
drop ± for brevity. The equation of motion 


v ( T o) = vfo) = i^/cv (to) (87) 

V w ( T o) 
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and hence 


H(to) 



w(t 0 ) 


( 88 ) 


where lo(t 0 ) = 12v / 2na 2 (r 0 ). 

Now, if we assume |0) to be the initial quan¬ 
tum state, and take into account the both 
modes ‘±’, then the number of particles N at 
some later time t± can be found with the for¬ 
mula: 


<0|iV(n) +1|0) 


<0|H + (n) + H-(n)|0> 


w(n) 


(89) 


|^| 2 (T 1 )+288n 2 a 4 (r 1 )|n| 2 (r 1 ) 

24V2na 2 (n) 


What remains to be solved is (83) with initial 
data (87). 




C. Numerical results 

We fix units so that Ai = 1 = A 2 and also 
n = 1. So the universe contracts from a = 00 , 
bounces at a = 1 and re-expands to a = 00 . We 
set the initial data as a(r 0 ) = 10 4 . In Fig. (4) we 
plot the number of particles and the scale factor 
versus time. We also find that varying initial a 
from 10 4 to 10 does not affect the number of 
particles produced, which never exceeds 0.35. 


VI. CONCLUSIONS 

In the present article we have examined the 
quantum dynamics of the vacuum Bianchi IX 
model, the Mixmaster universe. The Hamilto¬ 
nian constraint, which consists of isotropic and 
anisotropic variables, has been quantized. This 
split of variables is crucial both for implement¬ 
ing our procedure and interpreting the result. 
Suitable coherent states, namely the ACS, have 
been employed to obtain some insight into the 


FIG. 4. The evolution of the number of particles N 
and the scale factor a from the contracting phase 
through the semiclassic.al bounce to the expanding 
phase. We fix Ai = 1 = A 2 , n = 1, and o(r 0 ) = 10 4 . 

involved quantum dynamics of isotropic back¬ 
ground. Making use of adiabatic approxima¬ 
tion, we have identified the eigenstates for the 
oscillating anisotropy at its lowest excitation 
levels. Our procedure, developed by qualita¬ 
tive arguments and based on reasonable and 
tractable approximations, i.e. the replacement 
of the anisotropy potential with the harmonic 
potential and the Born-Oppenheimer-type evo¬ 
lution, is validated by considerations outside the 
adiabatic approximation. 

The main features of our quantum model 
are the following: (i) the singularity avoid¬ 
ance due to a repulsive term regularizing the 
singular spacetime geometry; (ii) the reduced 
contraction rate of the universe due to sup¬ 
pressed growth of the energy of anisotropy at 
the quantum level as it becomes frozen in a fixed 
quantum state; (iii) the stability of quantum 
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Friedmann-like state in the quantum Bianchi 
IX model both in the contraction and expan¬ 
sion phase. We emphasize that the anisotropic 
oscillations have a non-zero ground energy level, 
which means that there is no true quantum 
FRW state. 

The resolution of the singularity is due to the 
repulsive potential generated by the ACS quan¬ 
tization as in [20] . We note that our approach to 
‘quantizing singularity’ is very natural. It was 
developed within an appealing probabilistic in¬ 
terpretation of quantization procedures and, in 
a sense, it extends what is usually meant by 
canonical quantization (see Appendix A). Our 
approach is universal in the sense that it re¬ 
moves the singularity from anisotropic models, 
recently shown to be true also for Bianchi I 
models in [36]. 

In our framework, due to the Born- 
Oppenheimer approximation, the anisotropy 
degrees of freedom are assumed to be in a quan¬ 
tum eigenstate. This property leads to the ra¬ 
diation energy density. In the near future, we 
will focus on developing our scheme to include 
transitions between different energy levels by 
more detailed computations. The Born-Huang- 
like approximation is just a refined version of 
the adiabatic approximation as shown in Sub¬ 
section (IIIB). Qualitatively, the Born-Huang- 
like approximation does not change the behav¬ 
ior of the semi-classical Hamiltonian obtained 
within the Born-Oppenhcimer-like approxima¬ 
tion, because the former generates a supplemen¬ 
tary positive term oc q~ 2 that only renormalizes 
the term already present (due to our affine CS 
quantization). On the other hand, if we apply 
only canonical quantization to the system, then 
no term in q~ 2 is present at the beginning in the 
Born-Oppenheimer-like approximation. In that 
case the Born-Huang-like approximation gener¬ 
ates a new term in the semi-classical behavior. 
Since this repulsive term is responsible of the 
resolution of the singularity, we can say that 
with or without ACS quantization, the resolu¬ 
tion of the singularity holds true in the case of 
the Bianchi IX model within the framework of 
Born-Huang approximation. 


We emphasize that our approach is com¬ 
pletely different from the one of Misner’s, which 
was based on the so called steep wall ap¬ 
proximation, discussed in Appendix (B). In 
the steep wall approximation, there is no ten¬ 
dency for quantum probability to be peaked in 
the minimum of potential, and therefore there 
would probably be no quantum suppression of 
anisotropy as it occurs for the real potential re¬ 
vealed in the harmonic approximation of our pa¬ 
per. 

We have shown that the wave remains in 
its lowest energy states during the quantum 
phase. Even beyond adiabatic approximation 
there is no significant excitation of the wave’s 
energy level. It is interpreted that the quantum 
FRW universe, unlike its classical version, is dy¬ 
namically stable with respect to small isotropy 
perturbation. Thus, supplementing the FRW 
Hamiltonian with the zero-point energy origi¬ 
nating from quantization of anisotropic degrees 
of freedom provides a quantum version of the 
Friedmann model which could be used for a 
study of the earliest Universe. 
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Appendix A: Affine coherent state 
quantization 

Coherent state quantization is a particular 
approach pertaining to what is named in [30] 
integral quantization. When a group action is 
involved in the construction, one can insist on 
covariance aspects of the method. A detailed 
presentation of the subject is given in [30] and 
in Chapt. 11 of [37]. In this appendix, we give a 
short compendium of this approach before par¬ 
ticularizing to the integral quantization issued 
from the affine group representation. 
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1. Covariant integral quantizations 

Lie group representations [38] offers a wide 
range of possibilities for implementing integral 
quantization(s). Let G be a Lie group with left 
Haar measure d p{g), and let g H > U (g) be a uni¬ 
tary irreducible representation (UIR) of G in a 
Hilbert space R. Consider a bounded operator 
M on R and suppose that the operator 

R ■■= [ M (g) d/i {g ), M {g) := U (g) MU* {g) , 
JG 

. (A 1) 

is defined in a weak sense. From the left invari¬ 
ance of dp{g) we have 

U {g 0 ) RU f {go) = / M {g 0 g) dp {g) = R, 

JG 

(A2) 

so R commutes with all operators U{g), g G G. 
Thus, from Schur’s Lemma, R = CmI with 

cm = f tr (p 0 M {g))dn{g) , (A3) 

JG 

where the unit trace positive operator po is cho¬ 
sen in order to make the integral converge. This 
family of operators provides the resolution of 
the identity on R. 

I M {g) dv {g) = /, dv {g) := . 

J G C M 

(A4) 

and the subsequent quantization of complex¬ 
valued functions (or distributions, if well- 
defined) on G 

f ^ A f = I M{g) f{g)dv{g ), (A5) 

JG 

This linear map, function i —> operator in R, is 
covariant in the sense that 

U{g)A f U\g) = A U(g)f ■ (A6) 

In the case when / G L 2 {G , dp{g)), the quantity 
{ii{g)f){g') := f{g~ 1 g') is the regular represen¬ 
tation. 


A semi-classical analysis of the operator Af 
can be implemented through the study of lower 
symbols. Suppose that M is a density, i.e. non¬ 
negative unit-trace, operator M = p on R. 
Then the operators p{g) are also density, and 
this allows to build a new function f{g) as 

/(flO = A f ■= f tr{p{g) p{g’)) f{g')du{g '). 

JG 

(A? ) 

The map / i —> f is a generalization of the 
Berezin or heat kernel transform on G (see [39] 
and references therein). One observes that the 
value f{g) is the average about g of the original 
/ with respect to the probability distribution on 
(G, dv{g’)) defined by g' tr{p{g) p{g')). 

Let us consider the above procedure in the 
case of square integrablc UIR’s and rank one p. 
For a square-integrable UIR U for which \ip) is 
an admissible unit vector, i.e., 

c{ip) ■= I dp{g) | (^| U {g) |V>) | 2 < oo , (A8) 
Jg 

the resolution of the identity is obeyed by the 
coherent states \ip g ) = U{g) \if)), in a generalized 
sense, for the group G: 

[ p(g)dv {g) = I , dv {g) = , (A9) 

Jg c ( u ) 

p{g) = IV’s) {i> g I • 

2. The case of the affine group 

As the complex plane is viewed as the phase 
space for the motion of a particle on the line, 
the half-plane is viewed as the phase space for 
the motion of a particle on the half-line. Canon¬ 
ical quantization of the plane is covariant in the 
sense that it respects the translation symmetry 
of the plane through its Weyl-Heisenberg group 
extension. Our approach to the quantization of 
the half-plane follows the same principle of co- 
variance with respect the affine group structure 
of this geometry. 

Let the upper half-plane n + := {{q,p)\p € 
R, q > 0} be equipped with the measure dqdp. 
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Together with the multiplication 

(q,p)(qo,Po) = ( qqo,Po/q + p), q e K+,?eR, 

(Aio) 

the unity (1,0) and the inverse 

{q,p)~ l = (^,-qpj , (Ail) 

n+ is indeed viewed as the affine group Aff + (M) 
of the real line, and the measure dqdp is left- 
invariant with respect to this action. The affine 
group Aff + (R) has two non-equivalent UIR 
[40, 41]. Both are square integrable and this is 
the rationale behind continuous wavelet analy¬ 
sis (see references in [37]). The UIR U + = U is 
realized in the Hilbert space H = L 2 (M5_,dx): 

U(q,p)ip(x) = {e lpx /y/q)ip{x/q) . (A12) 

By adopting the integral quantization scheme 
described above, we restrict to the specific 
case of rank-one density operator or projec¬ 
tor p = \ip) (ip\ where is a unit-norm 
state in L 2 (ITj_, dx) which should be also in 
L 2 (Rlj_, dx/x) (also called “fiducial vector” or 
“wavelet”). The action of UIR U produces all 
affine coherent states, i.e. wavelets, defined as 
I q,p) =U(q,p)\ip). 

Due to the irreducibility and square- 

integrability of the UIR U, the corresponding 
quantization reads as 

f ^ Af = J n , 

. . (lA13) 

which arises from the resolution of the identity 

[ \q,p){q,p\ = i, (A14) 

J n+ 2ttc_i 

where 

r°° Hr 

<V= / \Hx)\ 2 -2+z ■ (A15) 

Jo x 

Thus, a necessary condition to have (A14) true 
is that c_ i < oo, which explains the needed 
square integrability of ip w.r.t. dx/x implies 


■0(0) = 0, a well-known requirement in wavelet 
analysis. 

The map (A13) is covariant with respect to 
the unitary affine action U : 

U(qo,Po)AfU ] {q 0 ,po) = A u( , 0jPo)/ , (A16) 

with 


(ii(q 0 ,Po)f)(q,p) = f ((qo,Po) 1 {q,p)) (A17) 

= /(A®( P -P0)) , 

it being the left regular representation of the 
affine group. In particular, this (fundamental) 
property is used to prove Eq. (51). 

To simplify, we pick a real fiducial vector. For 
the simplest functions, the affine CS quantiza¬ 
tion produces the following operators 

A p = -i-^- = p, A q p = C -^A qP j (A18) 
qf{x) = xf(x). 


Whereas A q = (co/c_i)q is self-adjoint, the op¬ 
erator p = A p is symmetric but has no self- 
adjoint extension. We check that this affine 
quantization is, up to a multiplicative constant, 
canonical, [A q , A p ] = Icq/c-iI. 

We obtain the exact canonical rule (i.e. A q = q 
and A p = p) by imposing cq = c_ i. This simply 
corresponds to a rescaling of the fiducial vector 
'll) as ipi{x) = tp(x / p) / ^fp, with p = Cq/c-\. In 
the remainder we assume that this rescaling has 
been done and therefore Co ('*/’) = c_ i (V 7 ) - 
The quantization of the product qp yields: 


A 


qp — 


c 0 qp + pq 
C- 1 2 


qp + pq 


= D. 


(A19) 


where D is the dilation generator. As one of 
the two generators (with q) of the UIR U of the 
affine group, it is essentially self-adjoint. 

The quantization of the kinetic energy gives 

A P 2 =p 2 + Kq ~ 2 , (A20) 

r oo l 

K = K(ip) = / ( i/j'(u)) 2 u -. 

Jo C_ 1 
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Therefore, wavelet quantization prevents a 
quantum free particle moving on the positive 
line from reaching the origin. It is well known 
that the operator p 2 = — d 2 /dx 2 in L 2 (R* + , dx) 
is not essentially self-adjoint, whereas the above 
regularized operator, defined on the domain of 
smooth function of compact support, is essen¬ 
tially self-adjoint for K >3/4 [42]. Then, quan¬ 
tum dynamics of the free motion is unique. 

As usual, the semi-classical aspects are in¬ 
cluded in the phase space. The quantum states 
and their dynamics have phase space represen¬ 
tations through wavelet symbols. For the state 
| (f>) one has 


®(q,p) = (g,p|</>)/\/27r, (A21) 

with the associated probability distribution on 
phase space given by 

pt(q,p) = ^yH<g,p| 0 >l 2 - (A22) 

Having the (energy) eigenstates of some quan¬ 
tum Hamiltonian H at our disposal, we can 
compute the time evolution 

p<t>(q,p,t) ■■= ( A23 ) 



For instance, any power of q is transformed into 
the same power up to a constant factor 


q p q? 


Cp-lC-p-2 p 
C-l 


(A27) 


Note that c _2 = 1 from the normalization of ip. 


We notice that q = coc_ 3 (c_i) _1 q. If we 
choose the fiducial vector such that c 0 = c_ 1 
in order to obtain the canonical rule, it remains 
q = c_ 3 <p Using a g-rescaling of the coherent 
states in the definition of the symbols Af like 
Af = (Xq,p\ Af |A q,p) allows to obtain q = q if 
we choose A = l/c_ 3 . 

Other important symbols are: 

P^P = P, (A28) 


2 . ~2 2 , C (l/) 

P I—>• p Z = p + 


2 ’ 


(A29) 


POO 

cWO = / W(x)f 

Jo 


Co 


1 -I- x ) da:. 

C-l 


for any state p. The map (A7) yielding lower 
symbols from classical / reads in the present 
case (supposing that Fubini’s theorem holds): 


f{q,p) 


1 

V^TTC-i 


r di 

Jo qq' 


e ip ^'- x ) Fp (q', X - X ') x 


iP 










where F p stands for the partial inverse Fourier 
transform 

F p (q,x) = -L f + °° e ipx f(q,p). (A25) 

V J —00 

For functions / depending on q only, expression 
(A24) simplifies to a lower symbol depending on 


qpi-> qp = C ° C 3 qp (A30) 

C-l 

Another interesting formula in the semi- 
classical context concerns the Fubini-Study 
metric derived from the symbol of total differen¬ 
tial d with respect to parameters q and p affine 
coherent states, 

pOO 

(q,p\d\q,p) = iqdp / (i>{x)) 2 xdx = iqdpc- 3 . 

Jo 

(A31) 

and from norm squared of d\q,p), 

do 2 

\\d\q,p)\\ 2 =c- 4q 2 dp 2 + L^, (A32) 

q 

r°° 1 

L = ^ dxx 2 {ip'{x)') 2 — - . 









2. Harmonic approximation 
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With Klauder’s notations [29] 


dcr 2 (q,p) :=2 [\\d\q,p}\\ 2 - \(q,p\d\q,p}\ 2 ] 

(A33) 

= 2 ^(c_ 4 - c 2 3 )g 2 dp 2 + L d ^ 2 j . 


Appendix B: Quantum anisotropies 


1. General setting 


The Hamiltonian H g of Eq. (21) reads 
H g = &2 P+ + 36^ 3 g 2/3 K(,g), (Bl) 


where 


K(/3) 



e 4 /3 + 


^ 2cosh(2\/3/3_) — e 6/3+ 




More explicitly, we have for H g the expression: 


Hg = 


2 A 2 n 2 ~ 

—r ? w). 


(B3) 


with 

£(<?) 


--A + X (q)v n (fd), 


A = d} + +dj_, 
(B4) 


x(<?) = 


I8A3 8/3 


V n (j3) possesses an absolute minimum for j3± = 
0, and near this minimum we have 


K(/3) = 8(/? 2 + p 2 _) + o(Pi). (B5) 


As mentioned above, V n (/3) and therefore £(q) 
possess the symmetry C 3 „. This group has 
three irreducible representations usually called 
Ai, A 2 and E. Therefore the eigenstates of £{q) 
can be classified according to these representa¬ 
tions. 


Using (B5) we obtain 

£{q)~- l -A + Z X {q){Pl + £-)■ (B6) 

Introducing the quantum numbers n± = 

0,1,..., corresponding to the independent har¬ 
monic Hamiltonians in /3 + and /3_, we deduce 
the harmonic approximation of the eigenvalues 
e(n+, n._) of £(q): 

e(?x+,n_) ~ 4 v / x(g)(?x + +n_ + 1), (B7) 

which gives the approximation of the eigenval¬ 
ues Ew(q) of H g , with N = n+ + n_, 

E N(q)~ + (B8) 


3. Steep wall approximation 


As mentioned above, taking into account the 
C 3 „ symmetry and the exponential walls of the 
potential, we can approximate V n (/3) by an equi¬ 
lateral triangular box as shown on the figure 3. 
The interest of this approximation is that it pre¬ 
serves the symmetry C 3 „ of the potential and it 
also possesses an explicit solution in terms of 
eigenstates and eigenvalues [43], [44], [45]. 

The size of the triangle is a free parameter 
that must be somehow adjusted, e.g., through 
some variational method. 

Let us denote by b the side length of the equi¬ 
lateral triangle box Tb of Fig. 5, and let us de¬ 
note by Ut the potential equal to 0 inside the 
triangle and equal to +oo outside. The sta¬ 
tionary Schrodinger equation — ip = e^ip 
with the Dirichlet boundary conditions has ex- 
plicit solution [43-45]. The eigenvalues ejn,n 
with m = 0, 1 , 2,... and n = 1 , 2 ,... are given 
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FIG. 5. Contour plot of V(/3) near its minimum. 
Red lines: the three C 3 ,, symmetry axis (3- = 0, 
f3+ = P-/V 3, f)+ = —/3-/y/3. A possible triangular 
box approximation in blue. 


by 2 


e 


(T) 

m,n 


87T 2 
3 b 2 



(B9) 


( T ) 

The ground state energy is e 0 {. The cor¬ 
responding normalized ground state wave func¬ 
tion 3 ^ 0,1 reads [43] 


^0,1 (£) 



2 sin 




47i ~/3+ , ^ \ , 
bV 3 3 / ' 

(BIO) 

cos—j, 


where /?± € T&. 

Assuming that the harmonic approximation 
(Sec. B 2) yields a fairly good approximation to 


the ground energy of £(q), the length parame¬ 
ter b can be estimated by imposing that the two 
ground states coincide up to some small correc¬ 
tion 


e 


(T) 

0,1 


Wx(q) ■ 


This relation yields 





(Bll) 


(B12) 


This leads to a new approximation e(m,n) of 
the eigenvalues of £ (q) as 


e(m, n) ~ 4 y/x(q) 



+ n + mn 


(B13) 


The eigenvalues Em of H g are still given for¬ 
mally by 


En( q) ~ ^nv/2^ 2 i?3 (N+l), (B14) 

but now N does not reduce to a simple integer. 
It is given by 

77l 2 

N= -|- n 2 + ran - 1, with (B15) 

3 

to = 0,1;, n = 1,2,... 

The last remark concerns the actual discrete 
spectrum of the Hamiltonian (Bl) (or (B4)). 
Clearly, it interpolates between the spectrum 
of the isotropic 2d-harmonic oscillator, which is 
linear in N and the spectrum of the triangular 
steep wall which behaves like squared integers. 
This is reminiscent of the ld-Schrodinger oper¬ 
ator with a symmetric Poschl-Teller potential 
for which the two types, harmonic and infinite 
square well, of approximations hold (see for in¬ 
stance [46]). In this article we have privileged 
the “harmonic” side of the anisotropy potential 
for obvious reasons, but it seems to us that a 
rigorous complete mathematical analysis of the 
spectral properties of the operator (Bl) is still 
lacking. 
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